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NOISE DUE TO ROTOR-TURBULENCE INTERACTION 
R. K. Amiet 

United Technologies Research Center- 


SUMMARY 


A procedure for calculating the noise due to turbulent inflow to a 
propeller or helicopter rotor in hover is summarized. The method is based on 
a calculation of noise produced by an airfoil moving in rectilinear motion 
through turbulence. At high frequency the predicted spectrum is broadband, 
while at low frequency the spectrum is peaked around multiples of blade passage 
frequency. This paper provides the results of a parametric study of the 
variation of the noise with rotor tip speed, blade number, chord, turbulence 
scale and directivity angle. A comparison of the theory with preliminary 
experimental measurements shows good agreement. 


INTRODUCTION 


One of. the potential sources of both rotor harmonic noise and broadband 
noise is that due to inflow turbulence. The turbulent velocity field produces 
a fluctuating angle of attack of the rotor blade leading to unsteady blade 
loading and the production of noise. If the axial length scale of the turbu- 
lence is such that a given eddy is chopped by more than one blade, the noise 
will tend to concentrate around multiples of blade passage frequency, while if 
a given eddy is chopped only once, there is no blade- to-blade correlation and 
the noise is broadband. 

The most rigorous method of treating this problem is to combine a model of 
the turbulence spectrum with airfoil response functions to determine the airfoil 
loading and an acoustical theory to determine the subsequent noise field. This 
is the approach used in references 1-3; the present study gives extensions and 
further results of the analysis presented in reference 1. This analysis is in 
turn an extension of the analysis given in reference 4 for an airfoil in 
rectilinear motion through turbulence which has been verified by extensive 
comparison with experiment in reference 5. 

The basis for the extension of the rectilinear motion theory to rotational 
motion is that so long as the acoustic frequency is somewhat greater than the 
rotational frequency of the rotor, the effects of rotation can be ignored and 
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the blade treated as being in rectilinear motion at each instant of time. The 
sound can then be calculated by averaging the calculated instantaneous spectrum 
for one revolution of the rotor. This is described more completely in referenc 
1. The present analysis treats the case of helicopter hover, vertical ascent o 
propeller forward flight with an axial velocity. An extension to the case of 
helicopter forward flight is expected to be completed in the near future. 

A preliminary comparison of experiment and theory shows excellent agree- 
ment between the two for the case of a model rotor with nonzero axial flow. 
Further comparisons between theory and experiment are presently being made. 

A recent review of this and other helicopter noise sources is given in 
the paper by George (ref. 6). 
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SYMBOLS 

Blade number 
Blade semichord 
Blade chord 
Sound speed 
Blade semispan 

Acoustic energy in a harmonic peak 
Energy spectrum of turbulence 
Acoustic frequency 
Airfoil response function 
Integrated airfoil response 

Chordwise turbulence wavenumber and spanwise turbulence wavenumber 

o>b/V 

Y 

Wavenumber range of energy-containing eddies 
Value of k for f = 1 Hz and = 0.628 
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Turbulence integral scale 


I Effective lift 
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Azimuthal, axial and relative blade Mach numbers 

Pressure jump on airfoil 

Effective radius of rotor 

Distance to observer 

Revolutions per second 

Acoustic spectrum 

Maximum spectrum level of a harmonic peak 

Rotor span 

Blade passage time 

Time 

Mean square turbulence velocity 
Azimuthal, axial and relative rotor velocities 
Unsteady velocity normal to rotor blade 
Vector position of observer 

Cartesian coordinate system fixed to rotor hub; nonrotating 
Cartesian coordinate system; y' along rotor span; x T along chord 
Azimuthal angle 
Gamma function 

* 

Angle of observer from rotor axis 

Acoustic wavelength 

Density 



Turbulence spectrum 
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0 ) Radian frequency of blade forces 

w 0 Doppler shifted frequency heard by observer 


THEORETICAL FORMULATION 
Notation and Summary of Analytical Method 

The geometry for the problem is given in figure 1 and is the same as that 
of reference 1. The rotor blade is divided into segments and the noise con- 
tribution of each segment is calculated. The relative velocity V r of a blade 
segment is related to the azimuthal velocity Vy and the axial velocity V z by 


V r 2 = V y 2 + V z 2 


(1) 


The basis for the analytical method is to utilize a noise calculation 
procedure for an airfoil in rectilinear motion through turbulence (see refs. 
4-5) . This calculation is applied to calculate the noise produced in the far- 
field by a segment of the rotor span as if that segment were in rectilinear 
motion. This is repeated for different azimuthal rotor positions, and an 
average spectrum is then obtained by integrating over the azimuth, remembering 
to account for the differing amounts of retarded time spent by the rotor in 
each azimuthal location. 


The final result for the far-field spectrum given in reference 1 is (note 
that the expression in reference 1 is for a single blade and is here multiplied 
by the number B of blades) 
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The function G represents the integrated airfoil response function and is 
given by 


wp 0 bz 
C 0 (r + M r x ') 2 



(3) 
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with 



g (£) e l ^ K y x/ y 1 d£ 


(A) 


Without the exponential factor in equation (5) , Jl would be just the airfoil 
lift per unit span. The exponential factor accounts for differences in 
propagation time to the observer from different points on the airfoil surface. 
Thus, Jt could be called an effective lift and it depends on observer location. 
A skewed gust of the form 


Wg= Wq 


i(cjt-Kx x'-Kyy') 

6 


(5) 


produces a pressure jump distribution AP on the airfoil. The airfoil response 
function g is then defined by the normalization 

g = APe ^ la,t Ky y V(27rP 0 v r W 0 ) (6) 

Analytical expressions for g can be found in references 7-8 and the evaluation 
Jl is discussed further in reference 1. 

In equation (2), T is the blade passage time, is the turbulence spectrum 

for the velocity component normal to the airfoil as written in the primed 
coordinate system, and K x and Ky are the chordwise and spanwise turbulence 
wavenumbers. The radian sound frequency to in an airfoil fixed frame is related 
to K x by 


K x = w/V r 


(7) 


The relevant Ky turbulence wavenumber is determined by the observer location; 
i.e.. 


Ky = K x M r y7(r + M r X') 


( 8 ) 


The radian frequency to 0 is the sound frequency measured in a nonrotating 
frame fixed to the rotor hub and so is related to to by a Doppler factor, i.e.. 


w _ I — M Z COS0 

w o | + M r x'/r 
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The airfoil semichord is b, d and s are the semispan and span, respectively, 
and R is the distance of the airfoil segment from the rotor axis. 


Energy in an Isolated Harmonic Peak 

When the harmonic narrow-band-random peaks in the spectrum are separated 
by deep troughs it is possible to derive a simple expression for the overall 
energy contained in each peak. To determine the energy in a peak, equation (2) 
must be integrated over w 0 between the troughs bordering the peak of interest. 
The peaks will be isolated (i.e., there will be a rapid drop off with w 0 on 
each side of a peak) if the X z or third argument of ^ varies rapidly enough 
with w 0 ; more explicitly, the peaks will be isolated if 


V r L» V 2 V r T 


( 10 ) 


where L is the turbulence integral scale. If this is the case, then the 
variation of K x and Ky with oo 0 can be neglected in the u) Q integral. The 
troughs in the sound spectrum are separated by 2 tt/T. Thus the range of inte- 
gration is to + 2ir/T where represents the trough to the left of the 
peak. Then noting that for any functions f (x) 
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the energy contained in a peak is found by integration of equation (2) to be 
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It will be noted that this is the same energy which would be calculated to be 
in the band < co < “]_ + 2 tt/T if there were no blade-to-blade correlation. 
Thus, blade-to-bla8e correlation acts to concentrate the acoustic energy 
around multiples of blade passage frequency, but it is not basically a method 
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of increasing the total acoustic energy radiated (unless one considers large 
correlation lengths and low frequencies such that blades are simultaneously 
correlated so that in the far field one adds pressures rather than intensities). 

For an observer on the centerline, a simple result can be obtained for 
the ratio of the sound energy in a discrete harmonic peak to the maximum peak 
level. For an observer on the propeller axis the y integration in equation (2) 
simply introduces the factor 2rr. At a blade passing harmonic, will be some 
multiple of 2 tt/T. Since the harmonic peaks are widely separated, then using 
equation (10) it can be seen that only the first term in the sum is needed. 

The maximum at a peak is then 


_ , 27 rBGV r _i . 

Spp max tu ° ,Z ^ = ~TVyV^ $ww(K x ,0,0) (13) 

In this equation, the parameter w to be used in G, and Ky is equal to m 0 
as expected since there is no Doppler shift for an observer on the centerline. 
This can be seen also from equation (9) by noting that 


x'= ^ [v y (x sinr-y cos/) - V z z] 


which becomes x' = -V z z/V r for an observer on the z axis. Likewise, for an 
observer on the z axis, equation (12) for the energy in a peak becomes 
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Taking the ratio of equations (13) and (14) gives 
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Introducing the Karman spectrum function 
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allows the Integral in equation (15) to be calculated giving 
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(17) 


This should be a useful result for calculating the overall intensity at a given 
directivity location. Strictly speaking it is only valid for an observer on 
the rotor axis, but it should provide meaningful estimates for other observer 
positions . 


CALCULATED RESULTS 


Several calculated spectra will be presented here for the purpose of 
illustrating the type of results expected, and for determining the dependence 
on the various parameters of the problem. In reference 1 it was shown that the 
results obtained by placing the entire span at an effective radius, R, equal to 
0.8 of the tip radius were not significantly different from the results obtained 
by integration over the blade span. Thus, the present calculations will use 
the effective radius approach. For the calculations presented here, unless 
specified otherwise, the following input parameters to the calculation will be 
assumed . 
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Figure 2 shows the effect of tip Mach number on the noise. At low 
frequencies, harmonics of blade passage frequency begin to stand out as shown 
by the dashed lines. In subsequent figures only the envelope of the peaks 
and troughs will be shown. The abscissa is a reduced frequency, Tc = wb/Vy , 
normalized by a reference value ko which is the value of Tc for f = 1 Hz and 
My = 0.628. Thus, since the smaller of the two tip velocity cases is half of 
the higher one, a given value of frequency for the smaller Vy case occurs at 
an abscissa which is twice the value for the large Vy case. The reason for 
normalizing the frequency abscissa in this manner is so that for a given 
abscissa value, both cases will show the results of an interaction with the 
same turbulence wavenumber component. In other words, for a given turbulence 
component, as the blade velocity is increased, the frequency of the sound 
produced will increase. 

Dipole noise is commonly associated with a velocity dependence. How- 
ever, this is based on certain assumptions which do not hold for the present 
problem, e.g., the turbulence intensity is here assumed to remain fixed as the 
rotor velocity increases, rather than increasing proportional to the rotor 
speed. In discussing the velocity dependence of the spectrum, the high fre- 
quency and the low frequency portion of the curve will be discussed separately. 

2 

For high frequency and for V r not too near 1, \Jl\ behaves as 1/V for 
fixed k. Thus, for large w and fixed k, G varies approximately as V r for an 
on-axis observer. For large w the summation in equation (2) can be replaced 
by an integral over k z , bringing the factor V z VyT/ (2irV r ) out front of the 
integral. The final result is that the high frequency portion of the spectrum 

O 

increases as V r at a given fixed Tc. The overall energy contained in the 
spectrum in the range from a given k value to infinity will increase as V r , 
however, since the plotted spectra represent the energy per unit Hertz rather 
than per unit Tc and there is a factor of 2 in the frequency ranges of the two 
curves. (Added note: equation (2) represents the energy per unit uj 0 , and a 

factor of 2 tt must be introduced to convert it to a per Hertz basis. The 
plotted spectra are on a per Hertz basis.) 

At low frequency \J.\ is independent of V r . Thus, for an on-axis observer 

O 

G varies as V r J for fixed k. The summation in equation (2) reduces to the 
single term 0) at a peak of a blade passing harmonic. For V z << V y , 

V r ^ Vy and the spectrum then is proportional to V r ^/T. The variation of Vy 

in figure 2 was obtained by varying the rotor rotational frequency, while that 

of figure 3 was obtained by varying the rotor effective radius at constant 

rotational frequency. Thus, T varies with Vy in figure 2 but not in figure 3; 

in figure 2 the amplitude of a harmonic peak varies as V r ^ and in figure 3 it 
varies as V r ^. In calculating the overall energy, it will be noted that halving 
Vy cuts the number of harmonic peaks in half. This together with the approx- 
imately 9 dB drop in the peaks gives approximately a 12 dB effect, the same as 
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for figure 2. Using equation (17) for the ratio of energy to harmonic peak 
amplitude shows that there is no additional V r dependence so that the final 
velocity dependence of the energy at low frequency is V r ^. 

These results are in accord with the results for an airfoil in recti- 
linear motion given by equation (17) of reference 9. If the turbulence 
intensity had been allowed to vary with rotor velocity, the velocity dependence 
of the acoustic energy would have been the more familiar results and at 
high and low frequency, respectively. 

Figures 4 and 5 show the effect of changing the number of rotor 
blades. At high frequency there is no blade-to-blade correlation and each 
rotor blade acts independently of the others. Thus, doubling B increases 
the noise by 3 dB . This is illustrated in figure 4 where 3 dB has been 
subtracted from the curve for B = 4 to illustrate that for high frequencies it 
then becomes coincident with the B = 2 curve. Figure 5 shows the effect on the 
envelope of the harmonics at low frequency. Note that the abscissa has a 
factor of B in it. Thus, for a given value of the abscissa, the two curves 
which have values of B differing by a factor of 2 will differ in frequency also 
by a factor of 2. Then every harmonic peak in one of the two plots will have 
a corresponding peak in the other plot. For this problem the first and sub- 
sequent blade passing harmonics are in the range k >> k e . Equations (16) then 
show that 'v- oo~H/3. Since the reduced frequency k is small, 2 is not a 
strong function of w and equation (3) shows that G ^ . The frequency of any 

given harmonic varies as B/T. Thus, from equation (13) Spp max ^ 

which is a weak dependence, and the levels of corresponding harmonics in figure 
5 would be expected to be about the same. This is seen to be the case. 

Equation (17) then shows that the energy contained in each of the harmonics 
for B = 4 will be twice that for B = 2. Thus, both the low and high frequency 
regimes give a 3 dB effect; that is, an overall doubling of the noise for a 
doubling of the blade number. It should be remembered, however, that the 
present calculation does not include any steady loading effects on the noise. 

If this were important, then decreasing the steady loading by increasing the 
blade number could have a significantly different effect. 

Figure 6 shows the effect of changes in blade chord on the noise. In 
equation (2) the chord appears explicitly only in the airfoil response function. 
When the sound wavelength becomes smaller than about 4 chords, the airfoil 
begins to act like an airfoil of infinite chord. Thus, changes in chord have 
little effect at high frequency as is noted in figure 6. For a given frequency, 
as the chord decreases the reduced frequency decreases to the point where 
eventually the airfoil is responding in a quasi-steady manner. In this regime 
the airfoil response is proportional to the chord, and the sound is proportional 
to the square of the chord. This can be seen by comparing the first harmonic 
of the C/10 curve (54 dB) with that of the C curve (72 dB, shown in figures 2-5); 
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i.e., a tenfold change in chord gives a 20 dB change in the noise at low 
frequency. Again it should be remembered that this ignores the effects on the 
other noise sources that changing the chord might have. 

Figure 7 shows the effect that variation in the turbulent length scale 
has. From equations (16) it will be noted that for an on-axis observer when 
k >> k e (10 f L >> Vr) ^ behaves as L“2/3. This gives a 6.67 dB decrease 
for each tenfold increase in L, a relation which is seen to hold in figure 7 
over most of the frequency range for all three curves. In fact, it is only 
below 100 Hz for the curve representing the smallest turbulence scale (L = C) 
that this relation does not hold. 

Figure 8 shows the effect of finite Mach number on directivity. The small 
Mach number case given has a directivity which is nearly that of a compact 
dipole. The directivity is not exactly that of a dipole since the rotor 
blades are aligned with the relative flow, rather than lying in the plane of 
rotation. All the curves are normalized to 0 dB at 0 = 0; i.e., on the 
upstream axis. The directivity plots for the higher Mach number cases tend to 
be pushed outward toward the rotor axis for most angles when compared with the 
low Mach number case. 

Figure 9 presents a comparison of theory with preliminary experimental 
results. The experiment was performed in the UTRC open-jet Acoustic Research 
Tunnel. The turbulence was generated by an upstream grid as in reference 5, 
and measurements verified that the turbulence was essentially isotropic. The 
experimental results are presented "as measured,” while the theoretical 
results include a correction to account for the presence of the tunnel open-jet 
shear layer through which the sound must pass before reaching the observer. 

It should be emphasized that there were no adjustable parameters in the theory 
which could be used to improve agreement between theory and experiment. The 
turbulence length scale and intensity were measured independently of the 
acoustic test. The large peak at blade-passage-frequency is due to the steady 
blade loading. Other than this steady loading effect which is not included in 
the present formulation, the agreement between theory and experiment appears 
accurate to within a few dB . Additional experiments are presently in progress 
to give further assessment of the theory. 


CONCLUDING REMARKS 


The change in the noise produced by a variation of the parameters 
affecting turbulence ingestion noise can be summarized as follows for an 
observer on the axis of the rotor: 
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(1) The acoustic energy has a dependence at low frequency and a 
dependence at high frequency. 

(2) The acoustic energy increases linearly with blade number at high frequency 
and nearly linearly at low frequency. 

(3) The effect of blade chord on noise is small for wavelengths greater than 
about 4 chords. For chords small compared to the wavelength, the noise is 
proportional to the square of the chord. 

(4) The noise spectrum varies with the -2/3 power of the turbulence integral 
scale except when L is comparable to or less than the chord and the 
frequency is low. 

(5) The directivity plots tend to be flattened toward the rotor plane when 
the Mach number is significant. 

The most effective method of reducing the noise appears to be through 
reduction of blade tip speed because of the rather strong velocity dependence 
of the noise, although not as strong as which one often expects for dipole 
noise. In lowering the tip velocity, other parameters would also have to be 
changed to maintain rotor thrust, and this could offset somewhat the beneficial 
effect of a decrease in tip speed. 

A comparison with preliminary experimental results gives strong support 
to the theory. A more comprehensive experimental assessment of theory is 
currently in progress. These experiments include the case of simulated forward 
flight in which the ambient flow is not axial. The present theory is being 
extended to treat this case. Also, although the calculations performed here 
were for the case of isotropic turbulence, the case of anisotropic turbulence 
can be treated if a satisfactory model of the turbulence spectrum is available. 
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OBSERVER 



Figure 1.- Geometry for rotating blade segment. 
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FAR-FIELD SPECTRUM (dB. RE 0.0002|j BAR) 











M = 0.628 

r 



0 -10 -20 -30 -20 -10 0 

FAR-FIELD SPECTRUM <dB, RE 0.0002|i BAR) 

Figure 8.— Effect of tip Mach number on directivity. 
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Figure 9.- Comparison of theory and experiment. M = 0.097; M 
c/L = 1.1842; B = 4, c Q /c « 6022.5; (u 2 /V = 0.067; R/c 
r/c = 4.3; r/c = 53.3; s/c = 5.33; RPS = 94.0; 9 = 124.85. 
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